function vA = vec(A)
vA = A(:);
end

